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ABSTRACT 


A Higher Harmonic Optimal Controller 
to 

Optimise Rotorcraft Aeromechanlcal Behaviour 


by 

Jane Anne Leyland, PhD/AMES 
NASA-Ames Research Centre 
Moffett Field, California 


Three methods to optimise rotorcraft aeromechanical behaviour for those cases where the rotorcraft 
plant can be adequately represented by a linear model with a system matrix, were identified and 
implemented in a 'stand-alone code. For convenience, the helicopter ^vibration ™Ju^n proWem wm 
selected as the subject problem for this investigation. These methods determine the optimal control 
vector which minimises the vibration metric subject to constraints at discrete time ^ints. These 
methods differ from the commonly used non-optimal constraint penalty methods such » those 
employed by conventional controllers (e.g., the Deterministic, Cautious, and Dual c ° ntrol '® r ^J" V'®{ 
the constraints are handled as actual constraints to an optimisation problem rather than i as just 
additional terms in the performance index. The first method is to use a Non-l'near 
alaorithm to solve the problem directly. The second method is to solve the full set of non-linear 
equations which define the necessary conditions for optimality. The third method is to solve each of 
the possible reduced sets of equations defining the necessary conditions for optima W 
constraints are pre-selected to be either active or inactive, and then to S| mply select the best 
solution. The effects of manoeuvres and aeroelasticity on the systems matrix are modelled by using 
a "proportional navigation" type pseudo-random pseudo-row-dependency scheme to define the 
systems matrix. 


Cases run to date indicate that the first method of solution (i.e., the direct optimisation of t he control 
vector subject to constraints, herein referred to as the "Optimal Controller") is reliable, robust, and, 
easiest to use. The algorithm employed for direct optimisation is particularly suitable to this problem 
since although it is designed to solve the general non-linear programming problem, it employs a 
successive quadratic programming method to solve this more general problem. Since method 
initially estimates the Hessian and then updates it successively as the quadratic soluhons are 
updated and since the only difference between the problem of interest and the standard Quadratic 
Programming Problem is that the constraints are quadratic rather than linear. < h e Hessian is invariam 
to tne optimisation process and is analytically known. Correspondingly, mod 'b^ a *'ons to this 
successive quadratic programming method which will enhance its reliability, overall robustness, and 
speed, appear feasible for this particular problem of interest. 


The second and third methods successfully solved systems of non-linear equations which define 
necessary conditions for optimality. It doesn't appear to be practical to use these methods by 
themselves however, because there exist many (perhaps infinite) solutions to these equations and 
no real way to recognise the solution yielding the global minimum or to even guarantee convergence 
to this solution if it were known. The second method can however, be used to verify ‘hat the 
necessary conditions are satisfied when the first method is employed. The use of the serond 
method for verification after the first method obtains a solution was made an option in the stand- 
alone code. 


The Deterministic Controller was added to the stand-alone eerie to Comrolle! 

romnarison Ootions to directly optimise the weighting coefficients of the Deterministic controller 
erthe P / in a spSSd r& to one aether (i.e., the "Conventional Controller" as referred to .herein) or 
individually (i.e., the "Optimised Conventional Controller" as referred to herein) whilst ^^^"19-^® 
constraints were provided as a means to obtain a more representative and meaningful comparison 
o° n <ronTrol?ers and to help access the relative merit of solving the actual optimisation i problem. 
Occasionally a mathematical conditioning problem occurs with the direct optimisation of the 
individual weighting coefficients of the Optimised Conventional Controller which causes a numerical 
overflow and subsequent error termination. This occurs when the optimisation process converges to 
a "solution" which has a harmonic phase angle with an ntt/2 value. In this case the v^ue ofone 
the associated harmonic coefficients approaches zero whilst the value of the associated weighting 
coefficient has no finite upper bound. 

Cases run to date indicate that the performance (i.e., the reduction of the vibration metric) of ‘he 
Ootimal Controller was superior to that of both the Conventional and the Optimised Conventional 
Sntro^ theory, the Optimal Controller yielded a ze ra vf rajort . metric <or 

sauare non-sinqular T-Matrices (i.e., when the number of measurements equals the number of 

controls with no 9 redundancy). As expected the. Opti m 'f " a?^n^?leTfo r ^l^les^xcept 
the Conventional Controller in performance, but inferior to the Optima Ctontrolle ’ a ^ x< l e P‘ 
a few degenerate cases where the performance of the three controllers was essentially equal. The 
performance gap was widest for square non-singular T-Matrices. 






Nomenclature 


A diagonal l.u.b. constraint limit matrix with dimension (' 2 Lx ^). The 
diagonal elements are the A n values in Eq (13). 

A coefficient of the scalar quadratic term in Relationship (34). 

a general vector with dimension (/xl). 

a coefficient of the scalar quadratic term in Eq (72). 

A n constraint limit for the n-th harmonic constraint equivalent to {l.u.b) y n in 
Relationship (13). 

A-, equality constraint linear coefficient matrix with dimension (M*M) in 
Relationship (40). 

A 2 inequality constraint linear coefficient matrix with dimension (M*M) in 
Relationship (41 ). 

B diagonal constraint weighting matrix in Relationship (31) with dimension 
(MxM). B is normally the identity matrix I M . 

B coefficient of the scalar linear term in Relationship (34). 

B general matrix with dimension (r x s). 

b general vector with dimension (^ x l). 

b coefficient of the scalar linear term in Eq (72). 

h equality constraint value vector with dimension (M*l) in Relationship (40). 

bi inequality constraint value vector with dimension (M*l) in Relationship 
(41). 

C proportionality constant e[0,l] for Eq (8) which defines the elements of 
the “flexible” part of the T-Matrix. 

C constant term in Relationship (34). 

C matrix of dimension (MxM) equal to F + F T . 

c scalar constant term in Eq (72). 

Qj (i,j)th element of C equal to fij + fji- 

D scalar l.u.b. in Relationship (35). 

F matrix of dimension (MxM) equal to T t W z T. 

f performance index value for the General Non-linear Programming 

Problem. 



Nomenclature (Continued) 


f(«) performance index function for the General Non-linear Programming 
Problem. 

fij (i,j)th element of F. 

ff quadratic coefficient matrix of dimension (M X M) in the performance index 
of the general QPP (Relationship (39)). 

G vector of dimension (M*l) equal to [zJw z t] • 

G T linear coefficient row vector of dimension (l*M) in the performance index 
of the general QPP (Relationship (39)). 

g.l.b. greatest lower bound. 

H scalar equal to Zjw z T 0 - 

M scalar in the performance index of the general QPP (Relationship (39)). 

HHC higher harmonic control. 

HVRP helicopter vibration reduction problem. 
i duty cycle number. 
i index for the elements of the 0-Vector. 

I M identity matrix of dimension (M x M). 

pseudo-identity tensor of rank three and dimension 

J performance index value. 

j(.) performance index function. 

J A augmented performance index for the Conventional Controller problem. 
J QPP performance index for the general QPP. 

J 1 performance index for the scalar non-linear programming problem. 

j 2 performance index for the scalar Conventional Controller problem. 

y augmented performance index value for the Min/Max calculus problem. 
^(•,*) augmented performance index function for the Min/Max calculus problem. 

k duty cycle number. 

L dimension of the “reference base set of sensors", the “base measurement 
dimension", the number of rows in the “core matrix”. 



Nomenclature (Continued) 

/ general dimension used for vectors and matrices. 

l.h.s. left hand side. 

l.u. b. least upper bound. 

M number of control variables to be optimised, dimension of the 0-Vector 
(the control vector), number of columns in the T-Matrix. 

N total number of sensors, dimension of the Z-Vector (the measurement 
vector), number of rows in the T-Matrix. 

n harmonic number. 

P number of equality constraints for the general non-linear programming 
problem. 

P general symmetric matrix of dimension (r x r) or (/ x /) . 

Q number of inequality constraints for the general non-linear programming 
problem. 

Q upper left sub-matrix of the partitioned Jacobian with dimension x M) . 

QPP quadratic programming problem. 

R upper right sub-matrix of the partitioned Jacobian with dimension x M) . 

R common value of harmonic magnitude constraint limits equal to ( l.u.b.) Yi 

for i = (n-1), n, (n + 1). 

Ry (i,j)th uniformly distributed pseudo-random number used to generate a 
T-Matrix. 

R k right-sub-matrix of the partitioned Jacobian for the k-th reduced set of 
equations with dimension (9 x 3). 

r general dimension used for vectors and matrices. 

f^ M Euclidean M-Space. 

Euclidean -Space. 

r.h.s. right hand side. 

s general dimension used for vectors and matrices. 

S lower left sub-matrix of the partitioned Jacobian with dimension (f x M) . 


iX 



Nomenclature (Continued) 


T quasi-static transfer matrix (the T-Matrix) as identified for the current duty 
cycle which is common to both Local and Global models, and which relates 

the Z- Vector to the 0-Vector with dimension (N*M). 

T general matrix with dimension (N*M). 

T lower right sub-matrix of the partitioned Jacobian with dimension 

quasi-static transfer matrix as identified for the current duty cycle for the 
Global model which relates the Z-Vector to the 0-Vector with dimension 
(N*M). 

T l quasi-static transfer matrix as identified for the current duty cycle for the 
Local model which relates the Z-Vector to the 0-Vector with dimension 

(N*M). 

Tij (i,j)th element of the T-Matrix. 
t time. 

VCS vibration control system. 

W z diagonal weighting coefficient matrix of the quadratic term in Z of 
dimension (N*N). 

W z general diagonal matrix of dimension (N*N) whose diagonal elements are 
all greater than zero. 

Wi i -th diagonal element of W z . 

Wke. diagonal weighting coefficient matrix of the quadratic term in A0 of 
dimension (M*M). 

W e diagonal weighting coefficient matrix of the quadratic term in 0 of 
dimension (M*M). 

We scalar weighting coefficient of the quadratic term in 0 in Relationship (36) 
for the scalar Conventional Controller problem. 

We, i-th diagonal element of W e . 

Y n amplitude of the n-th harmonic sinusoid. 

y n value at time t of the n-th harmonic sinusoid. 

Z measurement vector (the Z-Vector) with dimension (N*l). 

Z(«) measurement vector function with dimension (N x l). 


X 



Nomenclature (Continued) 


Zq actual measurement vector during the current duty cycle for the condition 
where no control is applied with dimension (N*l). 

Zq. i-th element of Zq . 

Z| i-th element of Z. 

Zj estimated measurement vector for the current duty cycle with dimension 
(N*l). 

Zj., actual measurement vector for the previous duty cycle with dimension 
(N*l). 

a slack variable vector with dimension (y * l). 

a k k-th element of a. 

a v diagonal slack variable matrix with dimension (f l y). 

r vector comprised of the a and X vectors with dimension (M*l) and equal 
to [a\X^ ■ 

i; vector with dimension (3*1) comprised of the k-th possible combination of 
elements from the a and X vectors as specified by Relationship (123). 

f 0 if i * j 

<%j Kronecker Delta = j ^ • 

Ad change in the value of 6 from that of the previous duty cycle with 

dimension (M*l) and equal to 0 k - 0 k -i for the k-th current duty cycle. 

A0 k increment in the k-th element of 6. 

£ general vector of dimension (M*l). 

T] general vector of dimension (N* 1) . 

T] 0 general vector of dimension (N x l). 

6 control vector (the 6 -Vector) with dimension (M x l). 

6 scalar control variable for the scalar non-linear programming and scalar 
Conventional Controller problems. 

6* solution value of 6. 


Nomenclature (Continued) 

ft coefficient of the n-th harmonic Cosine term. 

ft i-th element of 9. 

ft control vector (the ft -Vector) defined for the current or i-th duty cycle with 
dimension (M*l). 

ft.! control vector (the ft -Vector) defined for the previous or (i - 1) th duty cycle 
with dimension (M*l). 

ft^ l.u.b. constraint limit for 9 in the scalar non-linear programming problem. 

ft coefficient of the n-th harmonic Sine term. 

ft, constraint matrix form of 6 with dimension (Mxf). 

0 vector comprised of the 9 and a vectors with dimension * l) and 
equal to [fta] T . 

A adjoint vector comprised of % Lagrangian multipliers with dimension 

(fxi). 

A k k-th element of A . 

£ general vector of dimension (M*l) . 

£ (M*l) dimensional control vector counterpart of the 9 -Vector for the 

quadratic programming sub-problem of the successive quadratic 
programming method used to solve the general non-linear programming 
problem. 

2 vector comprised of Equations (1 03) and (1 00) with dimension * l) and 
equal to [Equation (103) i Equation (100)f which equals <f> c (Q) 

$(•) equality constraint vector function with dimension (P * l). 

</>(•) equality constraint matrix function with dimension (t x t)- 

0 (.) equality constraint vector function with dimension (t x1 ) formed by 
compression of the corresponding equality constraint matrix function 
$(•) into a column vector by the process described in Section 2.4.1 . 



k-th element of <p c (»). 



V> 

<D 

¥(•) 

¥(•) 

%(') 


%(•) 

V*) 

'F 
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Nomenclature (Continued) 

equality constraint vector function with dimension (M x l) for a QPP. 

vector comprised of the © and X vectors with dimension (2M*l) and 
equal to [©:A] T = [d:<x\X^ = [0T] T - 

vector comprised of the 6 and i; vectors with dimension (9 x l) and equal 

to [0irJ T - 

inequality constraint vector function with dimension (Q x l). 
inequality constraint matrix function with dimension (-f x -^)- 

inequality constraint vector function with dimension (-^ x l) formed by 
compression of the corresponding inequality constraint matrix function 
y/(») into a column vector by the process described in Section 2.4.1 . 

(i,j)th element of 

inequality constraint vector function with dimension (M*l) for a QPP. 

vector comprised of Equations (103), (100), and (104) with dimension 
(2M x l) and equal to [E • Equation (104)] which equals j • 

zero vector of dimension (M x l) . 
zero matrix of dimension 
zero vector of dimension (N x l). 


X i i i 


Nomenclature (Continued) 


Superscripts 

T matrix transposition 

* solution or unique solution to a problem. 

Subscripts 

C denotes a column matrix (vector) formed from a diagonal matrix by shifting 
the diagonal elements to a single column. 

G global model. 

i duty cycle or element position number, 

j element position number, 

k duty cycle or element position number. 

L local model. 

L dimension of the “reference base set of sensors”, the “base measurement 
dimension”, the number of rows in the ‘core matrix”. 

M dimension of the control vector (the 0-Vector). 

n harmonic number. 

N dimension of the measurement vector (the Z-Vector). 

QPP quadratic programming problem. 

0 value when no control is applied. 

1FD one-sided finite difference method. 

2FD two-sided finite difference method. 


1.0 INTRODUCTION 


Given the predicted growth in air transportation, the potential exists tor significant market 
niches for rotary wing subsonic vehicles. Technological advances which optimise rotorcraft 
aeromechanical behaviour can contribute significantly to their commercial development, 
acceptance, and sales. An example of the optimisation of rotorcraft aeromechanical 
behaviour which is of interest is the Helicopter Vibration Reduction Problem (HVRP). 
Although the HVRP was selected as the subject problem for this investigation, it is 
emphasised that the analysis described herein is applicable to all those rotorcraft 
aeromechanical behaviour optimisation problems for which the relationship between the 
harmonic control vector and the measurement vector can be adequately described by a linear 
model. 

The reduction of rotorcraft vibration and loads is an important means to extend the useful life 
of the vehicle and to improve its ride quality. Although vibration reduction can be 
accomplished by using passive dampers and/or tuned masses, active control has the 
potential to reduce vibration throughout a wider flight regime whilst requiring less additional 
weight to the aircraft than passive methods. 

Active control is achieved using a closed-loop "feedback” controller (Figures 1 and 2 and 
Reference 1 ). Typically, the measurement vector is defined from the measurements obtained 
with vibration sensors (e.g., accelerometers), the relationship between the measurement 
vector and the control vector (i.e., the helicopter plant matrix) is determined (i.e., plant matrix 
identification), and then a new control vector which will hopefully reduce the vibration metric is 
defined and used to control the system (e.g., rotor blade pitch). It is this last mentioned 
function, that is the selection of the new control vector, that is the subject of this research. 
Commonly used conventional classical controller schemes (e.g., the Deterministic, Cautious, 
and Dual Controllers described in References 1 and 2 ) employ a non-optimal constraint 
penalty method to define the new control vector, whilst the Optimal Controller which is the 
primary subject of this document employs a non-linear programming optimisation technique 
that accounts for the constraints in accordance with optimisation theory. 
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2.0 TECHNICAL 

The systems models used by the controllers described herein and their mathematical 
equivalency are presented first. Next, a heuristic "seat-of-the-pants" proportional navigation 
type technique used during controller performance verification to model aeroelastic effects for 
the Helicopter Vibration Reduction Problem (HVRP) is described. The equations used by 
both the conventional (Deterministic) and Optimal Controllers presented in Reference 1 are in 
general matrix/vector form. One of the purposes of this documentation is to present these 
equations in their expanded form which is amenable to computer coding. These equations 
together with their associated equations defining the necessary conditions for optimality and 
the required equations to solve the non-linear programming problem (i.e the analytic 
gradient, the Hessian, and the Jacobian) are then presented. The pseudo-random pseucto- 
row-dependency scheme used to generate the systems matrix (i.e., plant matrix) and the 
initial measurement vector to heuristically model the aeroelastic effects is describednext 
Then the Stand-Alone Optimal Controller System is described, and finally the results of a 
numerical study comparing the Optimal Controller with optimised conventional controllers for 
various "proportions" of aeroelasticity is presented. 
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2.1 General Controller Definition 

The general controller scheme (Figures 1 and 2) uses the measured state (i.e., the 
measurement vector) to define a new control vector which will hopefully reduce the v| bration 
metric during the next duty cycle. To accomplish this, it is necessary to have first assumed a 
plant model (i.e., a mathematical model relating the estimated measurement vector to the 
control vector). A linear model with a Systems Matrix (also commonly referred to as the Plant 
Matrix or the T-Matrix) that estimates the measurement vector is assumed for all the 
controllers mentioned herein. It is necessary to identify or define this systems matrix (T- 
Matrix identification) before the new control vector can be determined because the new 
control vector is, in general, dependent on this matrix. This identification of the systems 
matrix is in general, dependent on the actual measurement vector which defines the plant 
Sate of interest (i.e., not the estimated measurement vector defined by the assumed 
mathematical model) and is not a subject of this research. 
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2.1 .1 Systems Models of Controlled Vibration Response 

The mathematical model (i.e. , the plant model) relating the estimated measurement vector to 
the control vector is assumed to be linear for all the controllers mentioned herein. The Local 
and Global model and their mathematical equivalence are presented first, then the sensitivity 
of the plant matrix elements for higher harmonic Fourier models and their representation by 
uniformly pseudo-random numbers are discussed, and lastly the modelling of aeroelastic 
effects by means of a Proportional Navigation type scheme to introduce some degree of 
linear dependency to the rows of a non-square matrix is described. 
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2. 1 . 1 . 4 Sensitivity o# the T-Matrix Elements and Its Modelling Using Random Methods 

The helicopter applications of interest typically require the definition of the harmonic content 
of the responses of interest to the harmonics of a specified control. The responses are 
measured by various sensors (e.g., accelerometers) mounted on the helicopter and its rotor 
system. A Fast Fourier Transform (FFT) is then performed on this response data to define 
the coefficients of the associated Fourier Series which describe these responses of interest. 
The measurement (response) vector is typically comprised of the coefficients of selected 
higher harmonic sinusoidal terms of these Fourier Series. The specified control vector is 
typically comprised of selected mono-frequency Sine/Cosine coefficient tuples, or 
equivalently mono-frequency ampiitude/phase angle tuples, from those which describe the 
higher harmonic sinusoidal excitation of blade pitch. Once the measurement vector for the 
current duty cycle is defined, the T-Matrix can be identified by various methods (Reference 
3). The elements of fiie T-Matrix can, however, appear to have random values when 
compared with one another for a particular duty cycle and to appear to vary significantly 
between dtriy cycles because of the sens&vity of the components of the measurement vector 
(i.e., the selected higher harmonic coefficients of the Fourier Series which describe the 
responses of interest) to minor variations in the wave-form of the sensor measurements 
which include the effects of process and measurement noise. This apparent randomness 
appears to have uniform distribution rather than a distribution with cluster points or modes. 

This apparent randomness can be exploited to greatly simplify the T-Matrix definition process 
for controller feasibility and desirability analyses. The elements of the T-Matrix can be 
generated simply by using a uniform random distribution generator with a range of values 
appropriate for me problem rather than using a detailed rotor simulation for the purposes of 
the study reported herein. In general, the largest square sub-matrix of the T-Matrix wffl be 
non-singular (i.e., its rows and columns will be linearly independent) if its elements are 
generated in this manner. This fact provides a simple means to verify that the optimisation 
algorithm is working property for a square T-Matrix when the constraints are not active. In 
this case, the minimised performance index should be identically equal to zero (see Section 
2.3.2). 
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2.1 .1 .5 Modeling Aeroelastic Effects by Means of a "Proportional Navigation" Type Scheme 
to Introduce Some Degree of Linear Dependency to the Rows of a Non-Square T- 
Matrix 

In general, the elements of the control vector should be at least partially independent of one 

another. Fortunately this requirement is relatively easy to satisfy by the appropriate selection 

of the control variables to be optimised (e.g., mono-frequency Sine/Cosine coefficient tuples 
or ampitode/phase angle tuples from the expression defining the control). 

It is also desirable that the elements of the measurement vector be at least partially 
independent of one another to eliminate redundancy in the performance index and to reduce 
the computational burden resulting from a larger number of rows in the T-Matrix. FOr the 
HVRP the independence of the elements of measurement vector depends to some extent on 
the elasticity of the aircraft. For example, if a heScopter had purely translational motion and if 
two sets of three orthogonal accelerometers were firmly fixed in different locations in its 
fuselage, one would expect that the acceleration vector determined from each set of 
accelerometers would be the same if the airframe was perfectly ridged and there were no 
measurement errors. On the other hand if the airframe was not perfectly ridged (e.g., elastic), 
one would expect that the acceleration vector determined from each set of accelerometers 
would differ somewhat. In the first case, the two sets of accelerometers would be fuHy 
dependent while in the second case the two sets of accelerometers would exhibit some 
degree of independence from each other. 

The least number of possible independent measurements (i.e., the minimum number of 
measurement degrees-of-independence) occurs for an ideally perfect rigid aircraft whose 
motion is purely translational. For convenience, this least number of possible independent 
measurements for this "ideal aircraft" is referred to as the "base measurement dimension", 
and any reference set of sensors of "base measurement dimension" which spans the 
measurement degrees-of-independence space for this "ideal aircraft" is referred to as a 
"reference base set of sensors". Additional sensors mounted on this "ideal aircraft" will not 
provide any information in adtfition to that obtained from the "reference base set of sensors". 
Indeed, if any additional sensor is identically equivalent to, and has identically the same 
orientation as one of the reference base set sensors, the corresponding row of the T-Matrix 
for this redundant sensor would be identical to that correspondmg to its associated reference 
base set sensor even though the two sensors might be mounted at different locations on this 
"ideal aircraft". 

If the aircraft is flexible and "bends" under loads, sensors in addition to the "reference base 
set of sensors" will provide additional information, and the dimension of the measurement 
degrees-of-independence space will be greater than the "base measurement dimension . In 
general the greater the flexibifity of the aircraft, the greater the degree of independence 
additional measurements will have. This fact, together with the recognition that a perfectly 
rigid aircraft with pure translational motion (e.g., translational vibration) is the limiting case 
having the minimum number of independent measurements, suggests a simple means by 
which elastic and non-tra relational effects can be added to the simple random model defined 
in the previous section (i.e., Section 2.1 .1 .4). The degree of independence of the elastic and 
non-translational effects as measured by sensors in addition to the "reference base set of 
sensors’ can be simulated by a "proportional navigation" type scheme in which the required 
additional rows of the T-Matrix are defined as a linear combination of the row corresponding 
to one of the "reference base set sensors" and a new randomly generated row. 

The ideally perfect rigid aircraft whose motion is purely translational would be most efficiently 
modelled with a minimum size (L*M) T-Matrix where M is the minimum number of 
independent controls which comprise the Theta-Vector (i.e., the control vector) and L is the 
"base measurement dimension" (i.e., the least number of possible independent 
measurements for this “ideal aircraft"). For convenience this ideal minimum size T-Matnx is 
referred to as the "core matrix". Each row of this (L*M) "core matrix" corresponds to a 
specific sensor in the "reference base set of sensors". Each sensor in the "reference base set 
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of sensors* has a corresponding row of this "oore matrix", consequently this mapping is "one- 
to-one onto". Although it is usually the case that L is greater than or equal to M, it should be 
noted that the controller optimisation algorithm should be sufficiently robust to provide a 
usable non-catastropNc control update if sensors fail and L becomes less than M. An 
algorithm was devised to generate a synthetic T-Matrix which would simulate a flexible 
aircraft. This algorithm uses pseudo-random "uniformly distributed" numbers to define the 
rows of the more general ( N*M) T-Matrix system according to: 
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2.1.2 General Controller 


The Vibration Control System (VCS) described herein consists of a closed-loop controller 
(see Figures 1 and 2) which computes a new Higher Harmonic Control (HHC) vector to be 
used during the next duty cycle. This computation is designed to reduce vibration and is 
based on the latest measured state vector and the latest identified T-Matrix. The relationship 
between the measured state vector (i.e., the Z-Vector) and the control vector (i.e., the 0 - 
Vector is referred to as the "helicopter plant function" and can be highly non-linear. The 
•helicopter plant function" is that part of the actual helicopter itself which relates the measured 
state vector to the control vector. It is assumed that the “helicopter plant function" tor the 
controllers described herein can be approximated with sufficient accuracy between 
successive duty cycles by a linear relationship such as those described in Sections 2. 1.1.1 
and 2. 1 . 1 . 2 . It is emphasised that these Snear approximations define an "estimated" state 
vector based on a computed control vector rather than an actual measured state vector, and 
that in general the two state vectors will not be equivalent. The computation of a new control 
vector is, in general, dependent on the T-Matrix part of this linear relationship and 
consequently the T-Matrix must be known before this computation can be accomplished. A 
T-Matrix identification process such as one of those described in Reference 3 is usually 
employed to define the T-Matrix. It is assumed that the T-Matrix is identifiable and is 
consequently known for the subject research described herein. Identification processes 
themselves are not the subject of this research. The emphasis of this research is the 
computation of the new control vector for the next duty cycle given the T-Matrix. 
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22 . -Vbrsrtion Minimisation’ Via the Conventional Controller 

The Conventional Controller is simply a special case of the General Contr ^ i ^®® 

2 1.2 and Figure 2) in which a specific computational procedure is assumed for the 
determination of the control vector to be used during the next duty cycle. The assumed 
concutational scheme is similar to that of the Deterministic Controller which «s cfoscrtoedin 
References 1 and 2. The definition and disposal of constraints, the definition of foe 
performance index and its non-optimal augmentation with constraints, the 
problem to be solved and its necessary conditions for optimality, and the specific 
conputational procedure used to determine the new control vector are presented m Sections 

2.2.1 and 2.2.2. 



2.2.1 Conventional Controller Problem 


The Conventional Controller is the "classic" controller which has been in use for many years. 
Although sub-optimal, this controller has worked wen in many applications, particularly lor 
steady state operating systems whose control approaches a steady state value as the 
desired operating conditions are achieved. Instead of solving the "rear problem which seeks 
to minimise the "rear performance index subject to "rear constraints imposed on the system 
and the control vector, the Conventional Controller solves a simpler problem which does not 
in^ose any "rear constraints during minimisation. This is accomplished by adjoining a 
sp ec ffi c form of the constraints with weighting coefficients to the "rear performance index to 
form an augmented performance index, and then by directly minimising this augmented 
performance index. This constraint form is designed to be positive for all conditions and to 
approach zero from this positive "side" as the system approaches its desired steady state 
value. This technique of adjoining constraints to the performance index is referred to as 
"Internal Limiting" (References 1 and 2). The idea is that by minimising this augmented 
performance index the constraint functions will be minimised simultaneously with the “rear 
performance index. The motivation behind this approach is simply that a known analytic 
solution exists to this problem, while solution of the real problem has historically been much 
more difficult. The Conventional Controllers described in References 1 and 2 also employ 
what is referred to as "External Limiting" which is an after-the-fact imposition of constraints 
with corresponding adjustment of the control after the minimisation of the augmented 
performance index has been accomplished. The Conventional Controller used during this 
research and described herein does not employ this External Limiting. 

A convenient and simple form for the "rear part of the performance index which can provide 
an excellent measure of toe helicopter vixation is a simple quadratic metric of the individual 
vibration measurements; specifically; 
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It is enphasised that this choice of a performance index is but one of many possibilities. This 
performance index is relatively simple, representative of the control objective, and amenable 
to the pertinent mathematical derivations. 
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The consol vector lor the helicopter applications of interest is s^ed 

mono-frequency Sine/Cosine coefficient tuples from the sinusoidal terms wt^idescritie the 
hioher harmonic sinusoidal excitation of Wade pitch. The actual hardware fanftations rec*»re 

of the higher harmonic frequency sinusoids of interest (i.e., the rrwno-frequer^ Si^Cos^ 
coefficient tuples which form the control vector). For each frequency of interest, the 

corresponding sinusoid has the form: 

j.- 4 


where 
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The amplitude of the sinusoid Y n is then expressed (see Reference 4). 
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from which the "rear (or "actual") l.u.b. constraint can be written: 
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For convenience, this constraint is expressed in the form 
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where 
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Instead of attempting to solve the optimisation problem with the constraints expressed by 
relationship (13) treated as actual constraints to the problem, the Conventional Controller 
seeks to drive the harmonic coefficients in the control vector to zero by adding a weighted 
summation of their squares to the "real" performance index defined by equation (9) to form 
the augmented performance index to be minimised. This weighted summation is sometimes 
referred to as "internal limiting" and is expressed by 
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The Conventional Controllers described in References 1 and 2 also include a A 6 term of the 
form 
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This A 6 term is not included in the analysis reported herein for simplicity, and because its 
inclusion will not enhance the performance of the Conventional Controller. 

The Conventional Controller problem in its simple form is 

(<o 
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2.3 Vibration Minimisation as a Non-Linear Programming Problem 

The Optimal Controier is simply a special case of the General Controller (see Section 2.1 .2 
and Figure 2) in which a specific computational procedure is assumed for the determination 
of tire control vector to be used during the next duty cycle. The Helicopter Vibration 
Reduction Problem (HVRP) was posed as a non-linear programming problem and a 
su ccessi ve quad ratic programming method developed by Schittkowski, Stoer, and GiU et all 
(References 7 and 9 through 15) was employed to solve it. The definition of constraints, 
performance index, the problem to be solved, and the specific computational procedure used 
to determine the new control vector are presented in Sections 2.3.1 and 2.3.2. 
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2.3.1 General Non-Linear Programming Problem 


The general non-linear programming problem can be expressed in the form 
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The helicopter vibration reduction problem was described somewhat in the text prior to the 
definition the Conventional Controller Problem in Section 2.2.1. The performance index is 
that defined by Equation (9) and the control vector (i.e., the Theta-Vector) is comprised of 
mono-frequency Sine/Cosine coefficient tuples which typically describe the blade pitch 
sinusoids as per Equation (1 0). Since the control vector is comprised of tuples, its dimension 
is always even. It is assumed that there is an amplitude constraint of the form defined by 
Equation (13) for each mono-frequency tuple in the control vector, consequently the number 
of constraints is identically half the dimension of the control vector. The helicopter vibration 
reduction problem can be expressed as a non-linear programming problem in the form 
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/s rt-ferrcd H hereto 3-s //t “ Op/ tn*&l C*>m / roller- , 


11 




7£e_ yvu&sho* an'se^ as A U Ly LofLer fo fry bo So /i/e //e 
h*orz GonpItC&feJl fro* - 1 the^r proyra.*»mr*y problem Mf<*eA 
by ftebfroHsL/p fesj cjfe* fic dnve^b/'Oh^l Controller problem 
JUfi„ecf by &tUi**shf &) bas jprJuceJ ^Uf Us 

fo Lt s&-fcs/*cb>ry r^su/ls Itr A»^/ ye®rs . A yrzf/rc&/ 

Cxxvy&rrsoy, df a rca/sr exa^y/f of bofi, i^efLJs provide 

IvStCjlf ?»f° uLy ble Xo/u/so* fo Ht *tort Cot^pliC^'feJ 
yton - ft'r ear proyray*n*}*4 pro L/a*! cjouIcI b c pre. f-zr&h/? , 



Mtniinfs*'' J t - Ad +£d + C 
QeR" 

SUjecf fo: Q Z <D ^ f}^+\fi) - Q MaX 


uMerZ • 

A £ r «W D are rcs^/ar- c&^n^fs 

,y j j j 

A >■ o 


D > O 

s7y As 1 fie Ira./ performs ce )y>JtX 

Q ts fkc Ooftrol V&rl&h/t 

d HA< ts fit J.U.J, Co»ffr*t»f /*•■/ f*r Q 
M = 1 


( 3 ^) 
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fteJ&fton&L'/p %-of <3 S*C&/&r Cex^v&i-f fm*&l C*drv //*■ 

prololzi * » ?S ' 


Mini €. xj 2. ~ 


. J/ + 8 

{ * ^ , 04) 

I A $ +z?$-t-c-t- /i4 ^ 


-Su 4jec4 "/o . a/ ^ Cloias'! f^} w is pev~ 


uL 


ere 




/j fine I'/i'ferno'l limtfinij fc 

cons/r^'n $ 


J 2 / V 5 U * fpitJA-fcd p€r~f^y~Uh^Pic€ t y *d.e.\ 

td d /s Ur 5X:z,f&v- U}&gjL{ii«j CecJ^fClet , «/ ^or 7^e 

'Wef^aY I/tot-f/y* htrfa 

7Xe gr&pLic&l rzpncsettU'fioui o£ Yfcxe express ro*s for J{ y 

'Jzj <&nJl \Uq Q ts sbouh fa ff urt. Lc -fit Supers cri ft 

asferts k dr* oft fb vdut &{ s>o/u/t'oh of e tiler Op/i*,*! 

or Con\/e*t/ioti&l C&)4in)/ler problem ^ s.ppro prtd t _, rf c^*n 
cU&tr(y be seen &)*{ : 

% 

/. ~[tie So/ufion ptr jb\/yh a* c€ fsxJeK >J^ U Ur 

Godro//er problem fs /ess U±n Ue Solu'ftoi* perfor^^ucr 
i f)Jc-X d^ U fit dbh\Zen/iohd Cx>dro//eir prob Ufa , 


2 ? 



y A 

ji. %<L <?o C&drol V&r* *Lle @1 /© U* Of>h» &*Trt>H-er- 

proL/ejy, tS cflff-^r ^ / /A<5*m //e S’© fu f(0*n Gonial | feiridofe 
Q 4b //e Conventional CLnfroller ftolltJ** . 

3. /#<? 0 o/ufior, Control V&rt^iok 6 t<=> ft* C^fro/lvr 

ProUe-y* CZU* bit 0» fie Constraint th*if @/uja% I d *£*£££ 
aiol/e ifid'l* Ut SoLfion Cjonitrol I ravtaUt Q To ft* 
ConVenfieeai Grfroller ^ be U/»u> or aftoVC ftiS limtf 
O/Htx Only Cotn cider f ally on if 

4. 71 e akor e fire.*. aLserV*4ions ;«/'fs/ei ft *4 Here uill 
h aeeeraJ U discrepancies /« fte sJufton performance 
tnaleX StJ cedref. The V*lue of Jx Can never le /*■ ** 

JL*h -ft*/ of J( ^nJ. if /S possible, defending or Ut 

value of 4 J 8, C t an! kfa, ^d ao “'^ Le 

qredtr fU 6^ . ?„ f>r*dtc* , fie discrepancy in 

perform*** aU denial Can Le sttjn'fica nfy d^e. 
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2.Z.Z Solution 

Holtcaf^ef Vi'hrdto*. l^utucfron fro b lorn [UVR-P) defied 

by /Zel&(ionsL',y> ( 33 ) expressed to SfauJ&ra Form ts 

(37) 

(**) 


Mm, we J = d T FO + 1G t Q +U 

QslR 

SUjoxt to: Pfo) ’ A ~ &}& & f ^ On. 


U)l 


•err 


F f CiF^ <3n H ■^re die{-tn-e.d joy s (“LZ^J (23) J 

(2.4) J yespectt vety 
Q y [s defined. L>y (Z 2 ) 

/) a*i J, & ZtATt d*r(- ty, *d &CCorJi to fl&ir USZ-Cj C t *1 

F^u&d (<> * til) 

Om_ rs tbt letro iM&drix InitfL dt Hi ettS (“z.* 

2 

^(•) TS dt Co»s/r*;ni FiunC-ttOh M&frtx UlfL atfaenSto* 

(?<S) 


U rs useful convent e*»/ d/ ^*>/*,/ / 0 jDrrzset*d //Ve 

Simple , but h*por-fa*J , theorems dez.lt*q uidi, tuzHers o / 

Sy ^i«e^ positive <de£i* de *t ess p^r {ih e*i / fazdrtces • 

It /'$ sssu./ne^/ /-4a./ a-// ^a/r^ccJ vectors as used t* {de.se 

tii<£>v'ej>vis ^-<re c©^ (form a.4>/e &u.cl //a/ tLelr cJeh, ie*Is sure. 
fCa.( viuhojocrx . 


5 / 




atvcM 1 

w z ts 5 > d) a^o^i©/ of divnehs?** (Aj *A/) ^//d 

^ 3 ^ 0 <^d/ fi/fiwwA' &// ^ll!&/e(r //e*t ^ c^o 

j= /S Ve.c'&r -£ 

/y a*iJ /y &re [/ec/ors £ fj£^ 

~T ts & Vhdrt'x of divheyistoh (Aj X M) 

0 { d tS He Ojtro /cc/br oj. dt^cmstbh (fi) *t) 

^UJ As 6 Vc.C'foy- o£ dft*lCti<S foi* (Mxi) 


~~ft> 6 -h 


L yA 4 ^ > o ■€ $*\ ‘X * On 

or e^n/iAa/e**//^ ^4 As fosd/Ve Veft»tf? . 

*■ ftTUT]) >0 

^bc/ {/&l-ei*-fly [T l^f] /j fLsrrJjv-c 


h v %""%\\ut 0 0 ““°Ah 

o w z o — * & 

0 0 Ws • " ' ^ ft 


« o * 


o o o *■ ' ' j W(j)l 


32 




t T ^r = u 'tf* - • • • + m, 

y k > o Yk e {t j z i 3^ — W ] 

t T ^T >0 

or tJ 2 rs posrftn. De/^if? 

fk^i B 

dot* si (Air ^ ~ ~ '^L 

Utzn % -6 ^ 

Mk y M A VK { > O Yfc ^ W | £ 1 0„ 

u x = V + % f / * 

He. %■**£-*£*. rf 

fat = tyr*u [r$i - rfT T ]A/j t£] 

( = flTVzT] ) > 0 fj£ >j? M | ^ <? M 

or z^ui\/Jet~ily [T T ^zT] h foTijhe Dt/t»i/e Q,EX>. 
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Ptnce P ts syty/yie.i'rlc f P ~ p T 

“Then 

& T p& - b t p t s - [ps] t b - (e 

* # S P& syv* w e />-/' c # 


Pkeorenn 3 



firoo -f 

P fc Po^PiVC 0^1 in) /e h ~jI)€.On SA# / 

(i/g /S“ mmc'frtc S') nee Ulg Is Ji^c>h&l 

f~ ^5 ST^'/i'iAte/w'C hy Pteorenn F ~ F 

C - F+ F r - F-eF = £F 
» • C fs s 'yn*M*e-(rtc < 3 -* ^ C — - 2 _F 


34 



S)»ce. F~ rs /ost-fn/e Dtp*i-(e 

r T F f >0 \$* 0 M 

fc; = F\7.?]t, = u t fz, > o 

0 M 

• . £ rs l^esi-f/i/e Dtffatft Q.E.D. 


The-Or^ino 4~ 



Since P tS S'ytritH^rt c P ~ P^ 

z T pL ^[a T P]i - [a T P T ]l> = [pa ] T [ 

a T n = (nr*]j r = U T /°af 

S'^ce !jP& ts a sc&l&r-j L Pa - [& T Pa] 

a T PL, = ZYa 0 , £. p. 


3S 



S' 


7 ^ & depress /h<j yose~c( by $e.b&bi o*s?h( p (33) 

tf gtvo, -ftd 

Z($) ~ T9 +Z 0 (Sp*lto» (a.)) 


<^i <?/ 


T sytA&jr^ <£.^< 5 / non — s! h. tju 


No 


CjOHS 


ac'Z/Vf / i’isr! is 


Q r . B&, <■ A \f & £ So/u/tot. JJ»ac? 


It 


j-o fort-e* 

J a «••'*« Q 1 On ^~\J >0 f 0 +9' (M °M 


N 


frvob 

H>y Cra>to<Zr-s 7b<zaresi« / 

3 a u»t r ' &*H) m 9 - Z( 0 *) = T#*+Z a = 0 , 

Item J - [z(&*)^ z [z(d*)) = OjUzOu = (9 

iteorej^i i 

J >0 fZ£R u \?*Ou 

Or vzleh'f/y 

J >0 i 6£/z h \ e 1 e* 




/ieorc~M -5" provide r 3 Ch4 vyte&ns 7© Te-sY for error fa flc 

^Igorifinns &J code JestgutA h solve fU UVUP ptscA by fcddtonslly 
(33) , or eyujvalefdly by /feUd ton slips (37) &«A ( 2 $), If T ts non-Singylar 
SrUol square, and ff Ue Cetnsfrafals are no~f acfire J Hen He dgorifb*** 

sUJj yield d h*fa i muhn Ifoflct of \J &gutsl ~fo ^cW). 


fynresf l gallon df parlous wfloJs ~fo Solve Hit j/i/FP 3-s a toon- It ne ar 
progn* Vhyn fag problem* led To 4L< S*Jec4l on of fie highly SUCCESS fu,l 

Successive ^uaJra.'flc pr^gr^unyntng meflod developed, ly ScblHk ouiskij 
Sfoer y QH e/ a/ (ptfere nces 9 flrougf lf)d edited Solves die general 

lAondiueair programming fnbieim joy so Iving a sequence of related 
tfUad rail's. programming SU.b-p rob teens • Ont advan dage o f ibrs *»c iLd 
Is jb&rf yyadraflC programming problems Can be solved efficiently . 

A Very tmporfauf property of yuadraftc prvgra+nmtfag problems b H&t 

if fie yuodraf? C C&efftclenf fa&drlte fa fit* performance hndeX ts pbSiftVe 

Jitffaiftj fie pro l> jean has «£■ unigut Sofufto^ coblcl h y of course^ fire 

Cjl»bal s*>/ufrot* , 7>aj means fb&f & segment* of so/u/rons /© fir 

guadra&c programming Suj~ pnl/e ms C«)! I I Converge f~ iU f/oUj 

So/uft^j of fU general problem fa die //V/- 7t* m^tbt>4 rs 
CO^perfaf totally fap/emenfej fa ff? £%nJ-A/one Optimal (dm! roller 
£ysfem fSocfi'oh Xrf) by sdecfio*, of fie rtoSL Sysfen* (Pefe renC? ? ) 
of Codes (j.Z.j fie tMSL main driver voufthes DidCOh)F VhJCOA)& 
3 nJt -fbelr subroutines) (JLlcL U)ere designed /» s*o/*r yyiiie general 

hott-llhczr pngfeutiH'iKa finite using -fits h>eiit>4- 


3 ? 





vec /< 


o r~ 


Ts fie sea /ar t*deX 

f 5 ft\<L tjU'&Jira&i c Cocffrcle*! m\^4rtx iSrt 4 

dit^c-n j/on (Mxto) 

ts fit linear rou vector- L*i\tb 

dt*nt.*tS to” (i X M) 

/S' a scalar 

ts He e^uUity constraint linear Ce>z{{rcJe*'t 
yn&trtx u/ifl dime.* st «>/» (M 

faejusl'ity Constraint bnear £e>qfftcie*{ 
yn&tctx utrfL d)**en stort (MxM) 

ts He e/jud)ty Oons’tr&l d i/du-e. vector U>itl 

dimension (Mi l) 

is He in equably Const rat n! Value vector U'tHt 
dtn^ensi 0n (Mxl) 

Is Ux 2 ero \fe.c(or Ulifl di*it*sion (M* 0 
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tS Ht* Con fra / \ZCC~ft>Y~ UtH (Mx() 

ts Hz equality Oo»s~fr&t(*f £uncfto* t/ec/©r 
U>iH Aimers ran (Mxi) 

ts He C*>ns/r*ti*f vector 

CJiH cli'f*e»sfo* (Mx() 


755e ^u&Jr&irc £<aL- frol/emS UScA by He 

StxccessNe. ^u^Jr^hx p roy r&*n miry meHocf fo* 
by ust»<j C^pprOX! o f Hz ^Cr'^orm&tice. 

tmjex hhezr- approx s of Hz Comsfr^imf ^uref/ors 

a/ ea-c 4 sHp of He lier^Hon process, Ybese ^ fprox / ** dior s 
are obHineA lay Simple 0f He / <4 y &n<A fs 

JL»cHo»s /* fitU 4 to*s Lp (3o ) UtH Heir~ ^.pproprt ^Hfy 
jruMcddi z^a/r/X Tky/or Skates exps^st o^s, cjlerz tf fie 
^esj/a*, of £(+) ts i**f f °S i 7/Ve Jeftrift , He zJgorifL** 

^Jtfusls ){ Sb fL'J H /S'. He Suh- problem ts : 




tol\ 


erf 




V 




Is fit fosrfivt Jlefohl't-e &jppro of 

^ fie f/e.ss'/Ai^ of {(*) ev±U±'feJ a./ 0^ 

u)lfl\ dik*e+'Sfo r\ (j^)x.M) 

I ts fie Gradient °( H') cv*tuv!tJ a-/ Bjf- 
6^ I CSiJh o[i i**enSios\ (lx M) 

| ts fit <jj radical o {#■) tv&luaJeJ &•( 

Q k J fcijft J< t+ttnSi otn (FkM) 



A 


fs fie of pi*) evalira/ej 3/ Bjf 

UJl{L jlh*\enS}0 * (QkM) 


b ts He Gon'fr&l Vector £bcut/ erf*wf of () for 

fie Ou^dt^t'iC. prO^IT&l*h«i»J S'is.lo- pcolo l tin 

UlL dintnsto* (hxi) 

If opfi aS >?o/ £>o.LitveJ a/ //e C&nvpfefio'i of 

rf entire sfep j ffe Hess/a* rs upj&iej (&fere*ces 10 
//) ^ /ieA^ s/e^ to ^//e^/c^ 4 

/"/ /y o*ce nofed &*td &*tp>f&stsedf fl*d tf the tfp&dlrdtc. 

C&e.ff('ci& , '{ Wi ^ ikr* ^ T^ff perforins* ct t*JeX (f.t.j °3i th 

j&Uiomshp (3s) ffessta* ) ** & Z )) 

Is f°sjkv& fie* fleet Ts a unique So/ufh* (l.e.j Q 


<£or flthboAslips fej) j (to); (tl) j &*d ^ for lSi.lt ithsllfs 
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(^l-) f ( 4 %) j and fff)) to flit QPP Uoiicb t S j Off Course j i/e 

$UJ S'o futi Oh - 

Tie attraction a f a metiod based on Solving guadradiC 
programing suL-f>roUe*»s n fie c/osre resell lanoe flad 
“fLe. f(\/PP has ft> ffe QPP* Tbn^parrson 0 / "fie !d\/RP defined 

by iPelaftonslip^ ( 37 ) ( 3 $) tfiiL fit QfP daft* e d by 

Pt/ationsLys fejj zmd Oft) si outs /L-f: 

i. Tie performance. tndeX J of He IdVflP taj tded/Ca/fy 

He S&+*e ^u&dradic jx>\ri&, as Jqpp 0/ Pie QPP . 

X. The tdVHP bas ho eyuahdy Constraints j Or cyrutre>UnHy 
jif IS identically -fit / iert> l/c cfor A j is tdtndi cally fie 

^e.ro vna\r\x. * v». ^elaiJ&wvs^if p s (f°) ^^d C4 3) jor fie guadra~fic 

p ngramyning Snl-probtekn solved at cacb deration step* 

3 . Tie Only difference betwteu fie HV PP £nd a QPPfT fit 
Inequality Constraint function (jfrcL Cs quadratic for fit 
HVTP f/Pt laiiohsitp fa %)) &*d linear for fie Off (Pt/d'onslly 

m. 

i Since J 2 nd f'O) Of fit HVRP (pJdo«Mp ( 37 ) and ( 3 S)) 
are ijyadrbfics and k^onfh j fietr gradients are. known 
linear, and easy H eva/uade . If rs nded //a./ fie grad) cnf 

of a> ^uaArdic erasluded by fuo -side J finlft differ, 

is enact anal yd (rally (set SuL- Sccfion 2 . 5 \ 2 .i) &*d ben 
Use of gradienfs So evW*afe^ need[ nod detract t gradient 

accuracy. 


erenc-es 
enct 
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S.Si* J <!>(•) of fa fauf 3 fa fatlonslif S ( 3 ?) Z»Jl 

( 32 )) sure cs kh out* } fLctr Hcssra*s ave. 

i h v^-ri&h / &.*d k^oi^jn » According !y j tie ties stains i*e&d 
hot be ufddceA during fa process F* £uiurt 

Codes dedicated fa tfa Hl/fcP • 

(?, Since fat ana-Jratic Coe l(-taei*\ j*,afrrx F <?/ "/^e HVflP h 
positive Jeffaift / its Hess) a* ts dso post fat deft*! fa * 
->*J ea,cZ sr/e|0 0 / ^ }{erziton process yields tie Untrue 
So lu.{lon fa fhe CjjAoJid * c programming SuL- frol/ei*i snd 
fU UesSi&A he*-d not be adjusted fa be positive Jefondf . 

jilts U oFfyt Solution Ts He <^LfaI so tub to* Q^Ji He iterzflo* 
process Consequently Converges fa He global 'Solution 

in He limit. 

F^pericn cf fa d'de F* decodes fa{ fa Selected Success/i/^ 
^uodtefac jP<CQ cj r& t*j/h/'h a fae. tbod Ljorks Very Uell tuie^d *j 
fa is ire/ /fable ^ robust ^ r&kson&h/y c/ftctenfa &*d o.f?pe&cs 
to actually Converge to tfa ^loh^i Solution <?/ He 1{VP-P‘ 


jZ£ 

& Ti ts shown in Section CL.^t tfat the H&ssf&h &t ~f t 

equals fa+F T d.cL by Tteo/e*u 3 eyu<ds SLF tS 


positive deH»ifa 


del 



2.4 Vibration Minimisation as a Constrained Calculus Min/Max Problem 

Variations of the Optimal Controller which determine the control vector that satisfies the 
necessary conditions for optimality for the HVRP defined by Relationships (37) and (38) are 
conceptually possfole. It is emphasised however, that satisfaction of necessary conditions for 
optimality does not guarantee identification of the global solution. Optimal Controllers of this 
type (i e., controllers which attempt to solve the HVRP by determining a control vector which 
satisfies the necessary conditions for optimality) are simply special cases of the General 
Controller (see Section 2.1.2 and Figure 2) in which a specific computational procedureis 
assumed for the determination of the control vector to be used during the next duty cycle. 
The control vector solution to the non-linear equations which define the necessary conditions 
for optimality is determined using the Levenburg-Marquardt algorithm (References 6, 16, 17, 
18) The definition of calculus Min/Max HVRP problem and the necessary conditions for 
optimality are presented in Sections 2.4.1 and 2.4.2, the method of solution of the full set of 
non-linear equations which define the necessary conditions for optimality is described m 
Section 2.4.3., and the method of solution by selecting the best of solutions of each of the 
reduced sets of equations which define the necessary conditions for optimally when 
the constraints are preselected to be either active or inactive is described in Section 2.4.4. 
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< 2 ** 4 ,/ GtDKS’fm&h C&la^/us Mth/Max PnLU** 

It*. Hi/ fiP Si&td&d Fbrt* HephcJl Ly KeJ&fiousLips (%?) (3%) 

Cs ~hr&^^a>ci#)ex/ ~ft> %fahJ&ril Mt^/M&X C&lcu/uS forth by tdroducflDh 
of a f-uhcftoh of sl^clc V&rtdolt \Zccf&<r Qi f Q fb tneyw&lify 

Con.srfr&lh'i fuhcfio* fi-i uL^oL HsuS^in^s ti,fS ^ahcffo^ ~fc> 
d'tt esju&lify CoMsf^thf fuMcf'** <$(*). 7 ' s Accomf MU 
by #) to be 


<fi(e) = fte) -<£ol = o te) 

Uiik AiiWlSfo* ^ ) 

l^(*3 ft //<? i^e^uoIHy C&Asfr^'yii ^umcfton i^\&~(r'(X debited 

by PeU/iohs/^jo (z%j wfb dimehSi 0 * (^ X 

(X / v s Ue sbcle i^t rio,Lle vector u\iL Ji*e*s!o» fy*l) 
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Oly Cs sizzle \/2ft*Ue u&friK UlA Mimevsfr* 

govt&l eJetmeuts Qtrc fie frdii/tju&l 
s/aJr I r&n*Ues 0C kj § c^cd au e of 

tjbccL ts ossoci &{e~d ui fL a- a? tne^uaiiiy Oovsfr&fat 
i(')j k-i,^ 


, ^ 

■ • • -2. 




=■ 


(9 




0 


a 


a 
^ j 


Sub sit fu klA j fjzlSion'sb'if) (.3%) tv to fz^UoStOn ft^) ytdds fit 
desired \reSult' y specifically 

m - A -6p5^-<*^ - 0 (H) 

fact tU (or end nu kj k-tji,"- % t of 

#« ConS-fraJni tS ac^l/e (i . «•, [(9^2 0^ - /Ik); 

r/ lit aws/^J « M-i acfivt (f.e., </**.), the* Oi^fO. 

£e.ca. use of t^e j A, & , <W ^ aKie <$(G>) 

ts c/t&tjon^l. Mou{ ouy /oss fa tjeverdrfy y fits di^ovoi 
y»*htK $(&) 3 S defined Ly fyuo-fioh Off) ca*t ic* cai* pressed 

lurlo 3 column i+izfrix f Vector) (/?(©) Comprised of jud tit 


#. 7?e subscript C devotes a Column l/U&triX formed <(ro» 
a d't&flo^d matrix by fit's type of Co repression * 




ike. cl(a-fjoiA*( e/e^ci*fx (p(0j loy {be °(f~ di^cjOH&l 

^oro dew e>«7 s . ~]%e SVan^/a* rj[ Mt* /M&X C^(cu>lus fo&M of 

tie UVRP n 

MMmts c J = Q t F9 * 1G T 0 + // (*?) 

SuLjecf h i(&) -[A-^B£ r <0g t =O (4*) 
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Z.4-,1 Ajccesr&ry Co* dll to #s <(or Ophwxdify 

Tie heces's&ry Go^d'cho *s op~[u**di /y bz ve&djly 

J zrttfejl USi* g d&oji&rd Mfa/M&K CsJcm-Ius * Firsf t fo He 

p&r£o*~* l +^ MC 'f hndeX | [J(0), 4* c (©), M Mi* Mt ** 

//e e^uddy consir&td vendor ^(@) fo (Lt 

i/ul^X J(0} usl*g a*\ sdjoi'd \jedor \ ch**e*sto* ^.x/) 

£0 Lose c/ew cm/j fU \ fau.l'fsp ftejrs q&xM. o*e 

OjlytcL Corresponds io ^ Cofus'fra.tui'f . 


,<f.(s),x]. }[©>] -J(eh\ T 4>(0) to 




4 ? 



Ii rs U\d Since Uitse tjcyierd ho* - I their j 

more. {La* one So/ubon for & X c%n exi v j/ fc ~Ttes<$ ^oUr-f/'omai 

*SoUo»** Jo ho~f necessarily Correspond "/» reh/is/e j 

/c/ done fLe ^ jo hi l Minimum f of f^e HVR-f * » 

For He HVR-P } JU Gr&Jt^ p&rj of Gr*J[Q of |f®, X] 
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2.4.3 Solution of the Full Set of Equations Which Define the Necessary Conditions for 
Optimality 

Investigation of various methods to solve the non-linear simultanews s^em 
equations (i.e., Equations (50) and (51)) that define the necessary conditions ‘or^ir^for 
0 and A led to the selection of the IMSL system (Reference 6) of codes (i.e., th ® 
driver routines DNEQNF and DNEQNJ and their subroutines) which were designed to soive a 
quite general non-linear simultaneous system of algebrac equations by iSbrdi 
M arquardt version of the algorithm employed by MINPACI^tmes HYBRD and HVBRD1 
(References 16 and 17). These MINPACK routines are modifications of M- J- D. Powell's 
Hybrid method (Reference 18) which is in itself a variation of Newton s Method. 

Experience to date indicates that these routines work reasonably well for me HVRP i^hat 
thevdo reOabty converge to a "solution-. Unfortunately, in general there exist many (perhaps 
infinite) combinations of 0 and A which satisfy these equations which define the necessary 
conditions for optimality. Indeed, not only do solutions exist which correspond to relative 
extrema not the global minimum, but solutions can also exist which do not c °7 es P°^ t< ??J? | 
extremum at alt. The solution that the algorithm converges to is quite sensitive to the irututi 
(starting) estimate of 0 and A . There doesn't appear to be any practical way to select an 
iShestimate for © and A which will guarantee convergence to the desired global 
minimum unless of course, the solution is already known and is used as .J* 
estimate This method can, however, be used to verify that a previously determined solution 
obtained by solving the non-linear programming problem (Section 2 3.2) is u^eed a »W»n. 
Accordinaly options were provided in me stand-alone program (i.e. f options defined y 
. 11, 12, 13, and 14) which provide this verification to the solution obtained by solving the 

non-Bnear programming problem. 



2.4.4 Solution of a Reduced Set of Equations Which Define the Necessary Conditions for 
Optimality for Each Possible Combination of Active/Inactive Constraints 

Recognising that satisfaction of the full set equations (i.e., Equations (50) and (51)) that 
define the necessary conditions for optimality do not necessarily yield the global minimum, a 
more reliable method to obtain the global minimum was sought which was likewise based on 
the satisfaction of necessary conditions. The approach which was developed selects the 
"best" (i.e., that which yields the minimum) of the solutions of each possible set of reduced 
equations obtained when the constraints are preselected to be either active or inactive. Pre- 
selection of the activity of a constraint eliminates the corresponding switching equation (see 
Section 2.4.2) from the set of equations to be solved. Since there is a switching equation for 
each constraint, the number of switching equations is M/2 for the HVRP. Correspondingly, it 
is possible to reduce the dimension of the system of equations defining the necessary 
conditions for optimality from 2M for the full set of equations to 3M/2 for the reduced set by 
pre-selection of the activity for all the constraints. Since the number of arithmetic operations 
are proportional to the square of the dimension, an approximate reduction of 9/16 in 
computation time for solution of a reduced set of equations from that of a full set of equations 
is obtained. Unfortunately this method requires solution of several sets of equations, the 
amount being dependent on the M, and generally results in a net increase in computation 
time. But notwithstanding this increase in computation time, the solutions (accomplished 
using the same IMSL routines specified in Section 2.4.3) of lower dimension systems are 
more reliable in general. 
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Experience to date indicates that, as in the case of solution to the full set of equ at KDns, these 
marines work reasonably well for the HVRP in that they do reliably converge to a solution 
for each of the reduced sets of equations. But likewise, there also exists many (perhaps 
rti^eVcTmWn^s of © and A which satisfy these reduced sets of equations and not aU 
of them correspond to any extrema. Similarly, there doesnl appear to be any practical way to 
ceiArt an initial estimate for 0 and A which will guarantee convergence to the desired global 
minimum Since the computation times were greater for this technique than for solution of the 
a^nce mere did nSl appear to be any useful applied, on o th* 
technique s^h as providing verification that the necessary conditions were satisfied, further 
use of this method is not recommended. 
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2.5 


Expanded Equations Required for Solution of the Non-Linear Programming Problem 
and the Necessary Condtfons for Optimality 

The equations presented in the previous sections were, for the most part, in matrix/vector 
form and a bit more general than convenient for coding purposes. The corresponding 
expanded and detailed forms of these equations and those equations required for the 
computational processes described herein are presented subsequently. Specifically, the 
expanded constraint and performance index equations are presented first, followed by the 
analytic gradient and Hessian, and finally the necessary conditions for optimality and the 
analytic Jacobian. The proof of the interesting and useful fact that the numerically derived 
gradient of a quadratic function obtained by the conventional two-sided finite difference 
method is mathematically equivalent to the analytically derived gradient of the same function 
is presented right after the analytic gradient equations. 
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2.6 The T-Matrix Generation System 


The stand-alone code URAND was written to define the elements of T and Zo (see 
Subsections 2.1. 1.2 and 2.1. 1.3). Pseudo-random numbers were generated from a uniform 
distribution (Reference 8) and employed in this code to define these elements in accordance 
wi#r the random concepts described in Subsections 2.1. 1.4 and 2.1. 1.5. A listing of 
definitions of the principal parameters in URAND including all the input parameters and a 
brief description of the principal routines comprising URAND is presented in Appendix A. 
Listings of the VAX/VMS Command File used to execute URAND and the URAND FORTRAN 
routines themselves are presented in Appendix B. Input to, and output from, a sample case 
is presented in Appendix C. 
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2.7 The Stand- AJone Optimal Controller System 


The stand-alone code OPTIM was written as part of this research to provide a means to 
experiment with, to compare, and to evaluate the controller schemes described herein. 
Numerous controller options (see definition of (OPT in Append* D) were provided to specify 
the optimal controllers defined according to the theory described in Sections 2.3 and 2.4 and 
the pseudo-optimal controllers defined according to the theory described in Section 2.2 and 
variations thereof. Several IMSL MATH/LIBRARY routines (see Append* D and References 
5, 6, and 7) were employed in this code to accomplish tire required computations and to 
provide verification of analytic gradients and Jacobians. A listing of definitions of the principal 
parameters in OPTIM including afl the input parameters and a brief description of the principal 
routines comprising OPTIM is presented in Appendix D. Listings of the VAX/VMS Command 
Fie used to execute OPTIM and the OPTIM FORTRAN routines themselves are presented in 
Appendix E. Input to, and output from, a sample case is presented in Appendix F. 

Early on in the development of this code during the initial verification phases before the T- 
Matrix Generation System was developed, T-Matrices were obtained from actual 
CAMRAD/JA (Reference 19) closed loop HHC simulations of the BO-105 and S-76 rotor 
systems. HHC problems of various dimensions (i.e., various T-Matrix dimensions) were 
examined. A typical and illustrative example of the behaviour of the Optimal Controller as 
described in Section 2.3, was a 6-vector HHC applied to the four bladed S-76 rotor system 
witii a 6-vector measurement. The Control Vector for this case is: 



Control Victor 
Rotating Frame ( i3? ) 
Bladf Pitch 


The. Measuretoetrf Vector ts 
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3.0 CONCLUSIONS 


Tlvee optimisation methods to solve the helicopter vibration control problem were identified 
swd implemented. These methods attempt to determine the optimal control vector which 
minimises the vforation metric subject to constraints at discrete time points. These methods 
difrer from the commonly used non-optimal constraint penalty methods such as those 
employed by conventional controllers (e.g., the Deterministic, Cautious, and Dual Controllers) 
rn that the constraints are handled as actual constraints to an optimisation problem rather 
than as just additional terms in the performance index. The first method is to use a Non-linear 
Programming algorithm to solve the problem directly. The second method is to solve the fofl 
set of non-finear equations which define the necessary conditions for optimality. The third 
mrthod is to solve each of fire possible reduced sets of equations which define the necessary 
condrtions for optimally when the constraints are preselected to be either active or inactive 
and then to select the best solution. 


Cases run to date indicate that the first method of solution (i.e., the direct optimisation of the 
control vector subject to constraints herein referred to as the "Optimal Controller) is reliable, 
robust, and easiest to use. The algorithm employed for direct optimisation is particularly 
suitable to this problem since although it is designed to solve the general non-linear 
programming problem, it employs a successive quadratic programming method to solve this 
more general problem. Since this method initially estimates the Hessian of the performance 
mdex and constraint functions and then updates it successively as the quadratic 
programming solutions are updated, and since the only difference between the problem of 
interest and the standard Quadratic Programming Problem is that the constraints are 
quadratic rather than linear, the Hessian is invariant to the optimisation process and is 
anaiyticaftyknown. Correspondingly, modifications to this successive quadratic programming 
method which wiM enhance its reliability, overall robustness, and speed, appear feasible for 
this particular problem of interest. 

The second and third methods (i.e., solution of a full set or solution of reduced sets of non- 
linear equations which define the necessary conditions for optimality) successfully solved 
systems of non-finear equations defining the necessary conditions for optimality. It doesn't 
appear to be practical to use these methods by themselves, however, because there exist 
many (perhaps infinite) solutions to these equations and no real way to recognise the solution 
yielding the global minimum or to even guarantee convergence to this solution if it were 
known. The second method can, however, be used to verify that the necessary conditions 
are satisfied when the first method is employed. The use of the second method for 
verification after the first method obtains a solution was made an option in the stand-alone 
program. 

The conventional controller was investigated to provide a convenient means of comparison. 
In addition, options to directly optimise the weighting coefficients of the conventional 
controller either in a specified ratio to one another (i.e., the "Conventional Controller" as 
referred to herein) or individually (i.e., the "Optimised Conventional Controller as referred to 
herein) while satisfying the constraints was provided as a means to obtain a more 
representative and meaningful comparison of controllers and to help access the relative merit 
of solving the actual optimisation problem. Occasionally a mathematical conditioning problem 
occurs with the direct optimisation of the individual weighting coefficients of the conventional 
controller which causes a numerical overflow and subsequent error termination. This occurs 
when the "solution" has "zero" harmonic components whose associated weighting coefficients 
approach infinity as the optimisation process converges to the solution weighting coefficient 
vector. 

Casesrun to date indicate that the performance (i.e., the reduction of the vibration metric) of 
the Optimal Controller was superior to that of both the Conventional and the Optimised 
Conventional Controllers. In accordance with theory, the Optimal Controller yielded a zero 
vibration metric when no constraints are active for square non-singular T-Matrices (i.e., when 
the number of measurements equals the number of controls with no redundancy) or when the 


HZ 



proportionality constant (C in Equation (8) and input parameter RATIO to the T-Matrix 
G«%ration System) is 100 per cent. As expected, the Optimised Conventional Controfler was 
****** to Conventional Controller in performance, but inferior to the Optimal Controller 
for a I cases except a few degenerate cases where the performance of the three controllers 
was essentially equal. The performance gap was widest for those cases where the Optimal 
Controller yielded a zero vibration metric in accordance with theory. 
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Figure 1. Controlled vibration response. 



Figure 2. General scheme for the closed-loop 

vibration control system. 
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